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f^*) . Abstract 

We study the planar three-body problem and prove the absence of a complete set of complex mero- 
morphic first integrals in a neighborhood of the Lagrangian solution. 

^ . 1. Introduction 

l/~) \ The three-body problem is a mechanical system which consists of three mass points m 1; m^, m 3 which 

CN ■ attract each other according to the Newtonian law [16]. 

The practical importance of this problem arises from its applications to celestial mechanics: the bodies 
which constitute the solar system attract each other according to Newton's low, and the stability of this 

^**^ , system on a long period of time is a fundamental question. Although Sundman [21] gave a power series 

solution to the three-body problem in 1913, it was not useful in determining the growth of the system for 
long intervals of time. Chazy [3] proposed in 1922 the first general classification of motion as t — > oo. In 
view of the modern analysis [7] , this stability problem leads to the problem of integrability of a Hamiltonian 
system i.e. the existence of a full set of analytic first integrals in involution. Poincare [18] considered 
Hamiltonian functions H(z, ji) which in addition to z\, . . . , Z2n also depended analytically on a parameter fj, 
near \i = 0. His theorem states that under certain assumptions about H(z, 0), which are in general satisfied, 
^. ■ the Hamiltonian system corresponding to H(z,jjl) can have no integrals represented as convergent series in 

2n+ 1 variables z±, . . . , Z2n and n, other than the convergent series in H, \i. Based on this result he proved in 
1889 the non-integrability of the restricted three-body problem [22]. However, this theorem does not assert 
anything about a fixed parameter value \i. 

Bruns [2] showed in 1882 that the classical integrals are the only independent algebraic integrals of the 
problem of three bodies. His theorem has been extended by Painleve [17], who has shown that every integral 
of the problem of n bodies which involves the velocities algebraically (whether the coordinates are involved 
algebraically or not) is a combination of the classical integrals. 

However, citing [7] " One may agree with Winter [25] that these elegant negative results have no impor- 
tance in dynamics, since they do not take into account the peculiarities of the behavior of phase trajectories. 
As far as first integrals are concerned, locally, in a neighborhood of a non-singular point, a complete set of 
independent integrals always exists. Whether they are algebraic or transcendent depends explicitly on the 
choice of independent variables. Therefore, the problem of the existence of integrals makes sense only when 

rN ' it is considered in the whole phase space or in a neighborhood of the invariant set ..." 

Consider a complex-analytic symplectic manifold M, a holomorphic Hamiltonian vector field Xh on M 
and a non-equilibrium integral curve T C M . The nature of the relationship between the branching of 
solutions of a system of variational equations along T as functions of the complex time and the non-existence 
of first integrals of Xh goes back to the classical works of Kowalewskaya [6]. Ziglin [27] studied necessary 
conditions for an analytic Hamiltonian system with n > 1 degrees of freedom to possess n meromorphic 
independent first integrals in a sufficiently small neighborhood of the phase curve I\ One can consider the 
monodromy group G of the normal variational equations along T. The key idea was that n independent 
meromorphic integrals of Xh must induce n independent rational invariants for G. Then, in order that 
Hamilton's equations have the above first integrals, it is necessary that for any two non-resonant transforma- 
tions g, gi G G g must commute with gi. Although Ziglin formulated his result in terms of the monodromy 
group, it became quite recently [15,20] that much more could be achieved, under mild restrictions, by re- 
placing this with the differential Galois group. Namely, one should check if its identity component, under 
Zariski's topology, is abelian. 

The collinear three-body problem was proved to be non-integrable near triple collisions by Yoshida [26] 
based on Ziglin's analysis. 
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The present paper is devoted to the non-integrability of the planar three-body problem. 

In 1772 Lagrange [8] discovered the particular solution in which three bodies form an equilateral triangle 
and each body describes a conic. 

Moeckel [14] has shown that for a small angular momentum there exist orbits homoclinic to the Lagrangian 
elliptical orbits and heteroclinic between them. Consequently in this case the problem is not-integrable. 
Nevertheless, it was observed that for a large angular momentum and for certain masses of two bodies 
which are relatively small compared to the third one, the circular Lagrangian orbits are stable and, a priori, 
the system can be integrable near these solutions. Topan [23] found some examples of such transcendental 
integrals in certain configurations of the restricted three-body problem. 

Our approach consists of applying the methods related to [27,15] to the Lagrangian parabolic orbits. 
This means that we will study the integrability of the problem in a sufficiently small complex neighborhood 
of these solutions. 

The plan of the paper is follows. In Section 2, following Whittaker, we introduce the reductions of the 
planar three-body problem from the Hamiltonian system of 6 degrees of freedom to 3 degrees of freedom. 
Section 3 is devoted to a parametrization of the Lagrangian parabolic solution. 

Section 4 contains the normal variational equations along this solution. In Section 5 we study the monodromy 
group of these equations. In Section 6, applying the Ziglin's method, we prove that for the three-body problem 
there are no two additional meromorphic first integrals in a connected neighborhood of the Lagrangian 
parabolic solution (Theorems 6.2-6.3). Section 7 contains a dynamical interpretation of above theorems in 
connection with a theory of splitting and transverse intersection of asymptotic manifolds. 

2. The reduction of the problem 

Following Whittaker [24] let (xi, X 2 ) be the coordinates of mi, (x 3 , X4) the coordinates of m 2 , and (x§, x%) 

dx r 1 

the coordinates of 7713. Let y r = m/c— — , where k denotes the greatest integer in — (r + 1). The equations of 

motion are 

%-'£• t--& (r = 1 < 2 -- 6) ' 

where 

Hi = ^ O/i + vl) + 7. (vl + Vl) + o (Vb + % 2 ) - m 3 m 2 {(x 3 - x b f + (x 4 - x 6 ) 2 y l l 2 

Zmi 2m 2 2m3 

-77j 3 77Ji{(x5 - xi) 2 + (x 6 - x 2 ) 2 Y 1/2 - m 1 m 2 {{x 1 - x 3 ) 2 + (x 2 - xa) 2 }- 1 ' 2 . 

This is a Hamiltonian system with 6 degrees of freedom which admits 4 first integrals: 

T\ = Hi - the energy, 

T 2 = Vi + 2/3 + 2/5 j T 3 = y 2 + 2/4 + j/ 6 - the components of the impulse of the system, 

T4 = V\x 2 + 2/3X4 + J/5X6 — Xij/2 ~ x 32/4 ~ X5J/6 - the integral of angular momentum of the system. 

The system (2.1) can be transformed to a system with 4 degrees of freedom by the following canonical 
change (Poincare, 1896) 

_aw L am 

X r - ^ , 9r- g^ , (r-l,2,...,b), 

where 

(2.2) Wi = yih + y 2 l 2 + y 3 k + VaU + (2/1 + 2/3 + ys)h + (2/2 + 2/4 + Va)k- 

Here (l±, l 2 ) are the coordinates of mi relative to axes through m 3 parallel to the fixed axes, (I3, I4) are the 
coordinates of m 2 relative to the same axes, (l5,le) are the coordinates of m 3 relative to the original axes, 
(171,(72) are the components of impulse of mi, ((73,(74) are the components of impulse of m 2 , and (g5,ge) 
are the components of impulse of the system. It can be shown that in the system of the center of masses 
the corresponding equations for Z5, 1$, g$, ge disappear from the system and the reduced system takes the 
following form 

,o ^ dl r dH 2 dg r dH 2 



with the Hamiltonian 

rj M W 2 , 2n , M 2 2 , 2n , ! , , x ™3™2 TO1TO3 . mi™2 

#2 = -^-(fl-i + g 2 ) + — (# 3 + 9a) + (9193 + 5254 + , 

2 2 m 3 p 1 p 2 p 3 

where 

Pl = \J%+ll 92 = \jl\ + %, P3 = V(h-k) 2 + (l2-h) 2 , 

are the mutual distances of the bodies and M 4 = m^ + m^ , M 2 = m 2 + m^ . 

This system admits two first integrals in involution 
K\ = H 2 the energy, 
K 2 = g 2 h + g 4 l 3 + (?6^5 — gih — g 3 h — 55^6 = k - the integral of angular momentum. 

Let us suppose that the Hamiltonian system (2.3) possesses a first integral K different from K\ 2 . 

Definition 2.1 The first integral K of the system (2.3) is called meromorphic if it is rcprescntablc as a ratio 

K _ m,g) 
Q{i,gY 

where R, Q are analytic functions of the variables U, gi, 1 < i < 4. 

It can be shown [24] that the system (2.3) possesses an ignorable coordinate which will make possible a 
further reduction. 

Let us make the following canonical transformation 

(2.4) t-gJ, „--£, (, = 1,2,3,4), 

where 

W 2 = giqicosq 4 + g 2 q 4 smq 4 + g 3 (q 2 cosq 4 - q 3 sinq 4 ) + g4(q 2 smq 4 + q 3 cosq 4 ). 

Here qi is the distance m 3 nii; q 2 and q 3 are the projections of m 2 m 3 on, and perpendicular to mim 3 ; pi is 
the component of momentum of mi along m 3 m,i; p 2 and p 3 are the components of momentum of m 2 parallel 
and perpendicular to m 3 mi. 

One can write the new equations as follows 

dH 



(2.5) 


^i£ r \jj.± wy r 

dt dp r ' dt 


dq r ' 


(r = 1,2,3), 






and 












(2.5. a) 


dq 4 dH 
dt dp 4 ' 


dp 4 
dt 


= 0, 






with the Hamiltonian 












H -^H p ') 


M 2 , o 9n 1 


\ P1P2 


P3 p \ mim 3 

<?i J n 


m 3 m 2 

T2 


mim 2 
r 3 



P = P3(?2 -P2q 3 -Pi, 

where 

ri=qi, r 2 = yq 2 +q 3 , r 3 = yj (qi -q 2 ) 2 + q 2 , 

are the mutual distances of the bodies. 

Since p 4 = k = const the system (2.5) is a closed Hamiltonian system with 3 degrees of freedom. If this 
system is integrated then q 4 can be found by a quadrature from (2. 5. a). 

Proposition 2.2 If the Hamiltonian system (2.3) admits the full set of functionally independent meromorphic 
first integrals in involution {K±, K 2 , K 3 , K 4 } then the system (2.5) possesses two functionally independent 
additional first integrals {Hi, H 2 } which are meromorphic functions of the variables qi, pi, 1 < i < 3. 
This is the obvious consequence of the canonical change (2.4). 



3. A parametrization of the parabolic Lagrangian solution 

The equations (2.1) admit an exact solution discovered by Lagrange [8] in which the triangle formed by 
the three bodies is equilateral and the trajectories of the bodies are similar conies with one focus at the 
common barycenter. For the reduced form (2.5) the equality of the mutual distances gives 

(3-1) qi = q, q 2 = -, q 3 = — , 

where q = q(t) is an unknown function. Substituting (3.1) into (2.5) one can show that 

(3.2) Pi=P, p 2 =Ap^ , p 3 = Cp-\ , 

q q 

with p — p(t) unknown and A, B, C , D are the following constants 

m 2 (m 3 -m\) \/3kSim 2 m 3 \/3m 2 (mi+m 3 ) 

A — T, i £> = 77~ T, i ^ 



D 



miS 3 " S 2 S 3 mi S3 

km 2 (S 2 + m\m 2 — rrv£) 



S2S3 

where 

Si = mi + m 2 + m 3 , S 2 = mim 2 + m 2 m 3 + m 3 mi, S 3 = m 2 + 2m 3 . 

Substituting (3.1), (3.2) into the integral of energy H = h = const we obtain the following relation 
between q and p 

(3.3) ap 2 + -£ + - + — = h, 

q q q z 

where _ 

_ 2SiS 2 2 V / 3fcm 2 Si _ 2/c 2 Si(m| + m 2 m 3 + m|) 

a_ ^fs|' mi S| ' c --^ d - SIS~ 2 ' 

Moreover, from (2.5) we have 

, x da f A\ B 

(3.4) JL=( Ml + — )p+ 

at \ m 3 J m 3 q 

The equations (3.1), (3.2), (3.3), (3.4) define all Lagrangian particular solutions and contain two free param- 
eters: k and h. 

Consider the case of zero energy h = and k =/= 0. Then there exists a parabolic particular solution in 
the sense that the limit velocity goes to zero when the bodies approach infinity and each body describes a 
parabola. 

Putting w = pq one can find by using of (3.3) q, p as the functions of w 

(3.5) q = P(w), p 



where P(w) = —(aw 2 + bw + d)/c. 

Let M — C 6 be the complexified phase space of the system (2.5). Then (3.5), (3.1), (3.2) define a 
parametrized parabolic integral curve r G M with the parameter w € CP . 

4. The normal variational equations 

Let z = (qi, q 2 , q 3 ,pi,p 2 ,p 3 ), z e M. One can obtain the variational equations of the system (2.5) along 
the integral curve T 

(4.1) ^ = JH ZZ (T)C, CeT r M, 



where H zz is the Hessian matrix of Hamiltonian H at T and J is the 6x6 matrix 

-■(::)■ 

where E is the identity 3x3 matrix . 

These equations admit the linear first integral F = (£, H Z (T)), where H z = grad(H) and can be reduced 
on the normal 5-dimcnsional bundle G = T^M/TT of T . After the restriction of (4.1) on the surface F = 
we obtain normal variational equations (NVE) [27] which are the system of 4 equations 

(4.2) g =A(T) V , V eC\ 

where A is a 4 x 4 matrix depending on T. 

We can obtain NVE in the following natural way applying Whittaker's procedure [24] of reducing the 
order of the Hamiltonian system (2.5). 

Fixing the level of energy h — one can find p\ as a function of the other variables from the equation 
H(q,p) = which takes the following form 

a\p\ + hpi + ci = 0, 

where oi, b\, c\ are known functions depending onp2, P3, qi, 92, <?3- 
Solving this equation we get two solutions for p\ 



-fri + VA -&i-VA 
V\ = t, = K + and Pi = o = K -> 

where A = b\ ~ \a\C\. 

By substituting the Lagrangian solution given by (3.1), (3.2), (3.5) in these relations we choose the root 
Px = K- as corresponding to this solution. 

The functions q r {t), Pr(t), r — 2,3 satisfy the canonical equations 

(4 3) d<b_^dK_ dp L = _dK_ = ^ 

dq\ dp r ' dqi dq r ' 

where K = — K- and q\ is taken as the new time. 

The system (4.3) is a nonautonomous Hamiltonian system with 2 degrees of freedom which has the same 
integral curve T. Notice that K is not more a first integral. 

It is useful to pass now to the new time (q\ = q) —* w. From the formulas (3.3), (3.5) we have 

aw 2 + bw + d , 2aw + b , 

(4.4) q= , dq — dw. 

c c 

The resulting NVE (4.2) are obtained as the variational equations of the system (4.3) near the integral curve 
r and after the substitution (4.4) take the form 

(4.5) ^. = A(T) V , v eC\ 

where A is a 4 x 4 matrix whose elements are rational functions of w. 
We can represent A in the following block form 

i-C' T 

V -My 

where M\, M2, M3 are 2x2 matrices and Mj means the transposition of M; 




The matrix M\ is symmetric and has the following form 

M = 1 ( nn ni2 

1 S\L 2 Z 2 \ nis n 22 

where L(w) is the linear polynomial 

L = hw + h, 

and /i = 2S 2 , h = —V3mim 2 k. 

Z(w) is the following quadratic polynomial 

Z = Ziw 2 + z 2 w + z 3 , 

where zi = S 2 , z 2 = —V3mim 2 kS 2 , Z3 = k 2 m 2 (m 2 + m 2 rri3 + m 2 ). 
The coefficients mj have the form 

nn = Aiw 2 + A 2 w + A 3 , «i2 = A^w 2 + A 5 w + A e , n 22 — A 7 w 2 + A$w + A 9 , 

where Ai are constants depending on the masses mi, m 2 , 777,3 and fc. 
The matrix M 2 has the following expression 

^ ( X ° 

For the matrix M3 we have 



S 2 S\m\m 2 mz \ 1 



M 1 / ///,.., //',., 



miSlLZ \ 77721 ?7722 

where 

777ll = BlW 2 + -B277J + -B3, 77712 = B^VJ 2 + B5W + Bq, 77721 = B'jW 2 + -BgTi" + -B9 



777 2 2 = B m W 2 + BuW + B 



12- 



and Bj are constants depending on mi, 7772, 7773 and fc. 

The system (4.5) has four singular points w\, w 2 , W3, w^ in the complex plane: 



the infinity. 



the root of L 



u>i 



VO 777 1777 2 k 
W2 = ^^2—' 



(V37772 + iS 3 )kmi (V3777 2 - iS 3 )kmi 
(4 - 6) W3 = 2& ' W4 = 2& ' 

- the corresponding roots of the quadratic equation Z = where i 2 = — 1. 

Notice that the expressions for 7^2,3,4 have a rational form on the masses. 

The singularities Wi, 1 < i < 4 have a clear mechanical sense: toi corresponds to the motion of the bodies 
at infinity, w 2 defines the moment of the maximal approach. 

It is easy to see from (4.6) that if the angular momentum constant k — 0, then w 2 = W3 = W4, = and we 
have a triple collision of the bodies at the moment of time w — 0. If k =/= then by the lemma of Sundman 
there are no triple collisions in the real phase space and 77^4 become complex. 

Since the expression for p given in (3.5) becomes infinity when w — * 703 .4, formally, we can consider TU3 
and W4 as corresponding to the "complex" collisions which tend to 77; = as k — > 0. 

It was noted by Schaefke [19] that the equations (4.5) can be reduced to fuchsian form. 

In order to do it, consider the linear change of variables 

(4.7) 77 = Cx, 



where r\ = (771, 772, 773, t)a) t , x = (x\, X2, X3, X4) T and C — diag(LZ, LZ, 1, 1). 
In new variables the system (4.5) takes the following form 

(4.8) * = (J^ + J^ + .WL) X , xe <+, 

aw \ w — w 2 w — W3 w — w^J 

where A(k), B(k), C(k) are known constant 4x4 matrices depending on mi, m 2 , m 3 and k. 

Under the assumption k ^ we can exclude the parameter k from the system (4.8) by using the change 
of time w = kt. As a result, one obtains 

, An . dx ( A B C \ 

4.9 — = 1 1 x, 

y ' dt \t-to t-ti t-ti) 



where 



and 



V3mim 2 _ mi(^/3m 2 + iS 3 ) _ mi(^/3m 2 - iS 3 ) 

~ 2S 2 ' l ~ 2S 2 ' 2 ~ 2S 2 • 

A _ M(t ) B _ M(h) c= M{t 2 ) 



{to - ti)(*o - h) " (ti-t 2 )(ti-*b) («2-*i)(t2-to) 

Here, M(w) is the following matrix 

T dLZ 
3 

-L 2 Z 2 Mi 



MM= LZM >-lto E 




where one should put k = 1. 

The system (4.9) is defined on a connected Ricmann surface X = CP 1 /{to, ti, t 2 , 00}. 

It turns out that the matrix A is real and the matrices B = _R + i J, C = R — iJ are complex conjugate 
being i? and J real matrices. It will simplify matters further if we choose the units of masses as follows 

mi = en, m 2 = (3, my, = 1, < a < /3 < 1. 

In Appendix A we write the expressions for A, R, J with help of MAPLE. 

5. The monodromy group of the system (4.9) 

Let £(£) be a solution of the matrix equation (4.9) 

, N d ( A B C \ 
(5.1) -S= + + E, 

with the initial condition E(r) = I, t £ X where / is the unit 4x4 matrix. 

It can be continued along a closed path 7 with end points at r. We obtain the function E(i) which 
also satisfies (5.1). From linearity of (5.1) it follows that there exists a complex 4x4 matrix T 7 such that 
E(i) = E(i)T 7 . The set of matrices G — {T 7 } corresponding to all closed curves in A is a group. This 
group is called the monodromy group of the linear system (4.9). Let Ti be the elements of G corresponding 
to circuits around the singular points t = £,, i = 0, 1, 2. Then the monodromy group G is formed by To, Ti, 
T 2 . Denote by T^ G G the element corresponding to a circuit around the point t = 00. 

Lemma 5.1 T/ie following assertions about the monodromy group G hold 

a) Tq = I — is the unit matrix and 

(5-2) TiT 2 = T^ 1 . 



b) There exist two non-singular matrices U , V such that 



U- X T X U = V~ X T 2 V 



/ 1 1 \ 

10 

11 

\ o o o i I 



c) The matrix T^, has the following eigenvalues 

(5.3) Spectr(T 00 ) = {e 2 " Al , e 27TtX \ e - 2mXl , e" 2 ^ 2 }, 



wh 



c.r< 



(5.4) Xl = i + ^rs + Ve, \ 2 = l + ly/i3-V0, 



an 



d 



,S2 



= 144(l-3-|), S 1 = a + [3 + l, S 2 =a[3 + a + [3. 



Moreover, 



Spectr^)^}!,!,!,!}. 



Proof, a) The matrix A has the eigenvalues { — 1,-1,0,0}. Following the general theory of the linear 
differential equations let us write the general solution of the system (4.9) near the singular point t = to as 
follows 

x{t) = aXM + c 2 X 2 (t) + c 3 X 3 (t) + aX 4 (t), 

where C\ 4 G C are arbitrary constants and 

, , Xi(t) = -^=l-+ao + Oi(t-*b) + "- , X 2 (t) = — i- + bo + bi(t-t ) + --- , 

\°-°) t — to t — to 

X 3 (t) = co + ci(t - t ) + ■ ■ ■ , Xt(t) = do + di(t - t ) + ■ ■ ■ , 

where a,, 6,, c,, dj G C 4 are some constant vectors. 

By substituting (5.5) in (4.9) one can find a.j, 6,, Cj, rf, and show that the vectors Xl(£), X 2 (t), X 3 (t), 
X±(t) are functionally independent and meromorphic in a small neighborhood of the point t — to- This 
implies that the element To of the monodromy group G corresponding to a circuit around to is the unit 
matrix. Obviously we should have TqT\T 2 — T^ 1 . From this fact the relation (5.2) follows. 

b) The matrices B, C have the same eigenvalues {—2,-1,0,1}. It can be shown by a straightforward 
calculation that near the singular point t = t\ the general solution of the system (4.9) can be represented as 

x(t) = ciYi(t) + c 2 Y 2 {t) + c 3 Y 3 {t) + aYi(t), 

where ci ; ... ,4 G C are arbitrary constants and 

Y^t) = ei(t - ti) + e 2 (t - h) 2 + • • • , Y 2 =f Q + f 1 (t-t 1 ) + --- + Ciln(t - ti)Yi(t), 
Y 3 (t) = -^- + .90 + .91 (t - ti) + • • • , 

Y ^) = 77^2 + T^T + ■■■ + C Mt - ti)(/o + /i(t - ti) + ■ ■ ■ ) + C 3 ln(t - t^Y^t), 

[t-tiY t-ti 

where ej, /», g"j, hi G C are some constant vectors and C\, C 2 , C 3 are parameters depending on the masses 

a, (5. 



For C\, C 2 one can find 



d 



9/3a 3 (/3 + 2) 2 (a/3 + a + /3) 
4 (a + j3 + l) 3 ' 



C 2 = *Ci. 



The matrix S(i) = (Yi, Y^, ^3, Ki) represents the solution of the system (5.1) in a small neighborhood of 
the point t = t\. After going around of t\ we get £(£) = T,(t)M where 



M = 



I 1 2nid 27riC 3 \ 

10 

1 2niC 2 

\ 1 



Since C\ =/= 0, C2 9^ for a > 0, (3 > 0, there exists a non-singular matrix T such that 
(5.6) 



T~ X MT = 



/ 1 1 \ 

10 

11 

\ 1 I 



which is the Jordan form of M. 

The matrix 7i is similar to M and therefore has the same Jordan form (5.6). Repeating the analogous 
arguments for the matrix T 2 we deduce that the same assertion holds for the monodromy matrix T 2 . Notice 
that the existence of logarithmic branching near some Lagrangian solutions in three body problem was first 
observed by H. Block (1909) and J.F. Chazy (1918) ( see for instance [1]). 

c) Consider the matrix Aoo = — ( A + B + C) . Then there exists (see for example [4] ) a non-singular matrix 
W such that 

(5.7) ^ = W^e 2 ^ A ^W. 

Appendix A contains the expressions for the elements of the matrix Aoo . One can calculate its eigenvalues 

Spectr( J 4 00 ) = {Ai, A 2 ,3 - Ai,3 - A 2 }, 



where Ai j2 are given in (5.3). 
One can easy check that 

(5.8) 



o< Ve < 12, 



for all a > 0, (3 > 0. 

With the help of (5.7) we obtain for the eigenvalues of the matrix T^ the expression (5.3). 
Let us suppose now that Spectr(Too) = {1, 1, 1, 1}. Then according to (5.4) we obtain 



(5.9) 



13 + Ve 



Til, 



13- Ve-- 



n 2 , 



ni,n 2 € 1iZ. 



Hence, in view of (5.8), the number r = Vd is an integer < r < 11. The simple calculation shows that for 
these r the relations (5.9) arc not fulfilled. This implies that 

SpectrC^) 56 {1, 1, 1, 1}. 

The proof of Lemma 5.1 is completed. □ 



6. Nonexistence of additional meromorphic first integrals 

We call the planar three-body problem (2.1) meromorphically intcgrable near the Lagrangian parabolic 
solution r, defined in Section 3, if the corresponding Hamiltonian system (2.3) possesses a complete set of 
complex meromorphic first integrals (see Definition 2.1) in involution in a connected neighborhood of I\ 
Recall that equations (2.3) describe the motion of bodies in the system of the center of masses. 

From Proposition 2.2 it follows that in this case the system (2.5) admits two additional first integrals 
which are meromorphic and functionally independent in the same neighborhood. 

Theorem 6.1 For k ^ for the Hamiltonian system (2.5) there are no two functionally independent 
additional first integrals, meromorphic in a connected neighborhood of the Lagrangian parabolic solution T . 

Proof. Suppose that the Hamiltonian system (2.5) admits two functionally independent first integrals 
H i , H 2 , meromorphic in a connected neighborhood of the Lagrangian parabolic solution T and functionally 
independent together with H. According to Ziglin [27] in this case the NVE (4.5) have two functionally 
independent meromorphic integrals F\ , F2 which are single- valued in a complex neighborhood of the Ricmann 
surface T = CP 1 /{io,ii,i2, °o}- The linear system (4.9) was obtained from (4.5) by the linear change of 
variables (4.7) and the change of the time w = kt, k 7^ 0. Therefore, it possesses two functionally independent 
meromorphic integrals I\, I 2 . From this fact the following lemma is deduced 

Lemma 6.2 (Ziglin [27]) The monodromy group G of the system (4-9) has two rational, functionally 
independent invariants J\, J 2 . 

In appropriate coordinates, according to b) of Lemma 5.1, the monodromy transformation T\ can be 

written as follows 

/ 1 1 \ 

10 

11 

\ 1 J 

where / is the unit matrix and 



T, 



= I + D, 



(6.1) 



D = 



/ 1 \ 



1 

\ j 



For the monodromy matrix T2 one writes 



where 



T 2 = I + R, 



(6.2) 



R = VDV- 1 = 



( 01 a 2 a 3 a 4 ^ 

61 b 2 b 3 64 

c\ c 2 c 3 c 4 

\ d\ d 2 d 3 d 4 J 



with some unknowns a^, bi, Ci, di € C and a nonsingular matrix V. 
Let us input the following linear differential operators 

d d 

= ^2 7 VXi- , 

ox 1 0x3 



and 



A 



J2 a * 



dx\ 



5> 



dx 2 



E 



< j.i j 



dxv 



J2 d * 



dx 4 
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Lemma 6.3 Let J be a rational invariant of the monodromy group G, then the following relations hold 

8 J = 0, A J = 0. 



Proof. For an arbitrary ngNwe have T" = I + nD, hence J (T™x) = J (x + nDx) . Expanding the last 
expression in Taylor series we obtain 



(6.3) 



J {T?x) = J(x) + n8J(x) + J2 n l n{x), 



i=2 



where r^(x) are some rational functions. 

In view of J (T"x) = J(x) and the fact that J(x) is a rational function on x, the second term of (6.3) 
gives 8 J = 0. The relation A J = is deduced by analogy from the identity J (T2X) = J(x). □ 

Case (1). Assume that invariants J\ 1 J 2 depend on X2, x^ only. By Lemma 6.3 we have 



(6.4) 



A Jx = 0, A J 2 = 0. 



It can be verified that the equations (6.4) imply the conditions 6, = 0, di = 0, 1 < i < 4. Accordingly, 
the matrix R may be written 



(6.5) 



R 



I a x a 2 a 3 a 4 \ 



ci c 2 C 3 C4 

\ / 



One can find the characteristic polynomial P(X) = det(R — XI) of R 

(6.6) P(A) = A 4 - (01 + c 3 )A 3 + (01C3 - cia 3 )A 2 . 

In view of (6.1), (6.2) all eigenvalues of the matrix R are equal to 0, thus, with help of (6.6) we get 

(6.7) 01 + c 3 = 0, aic 3 = cia 3 . 
The matrix T\Ti takes the following form 



TiT 2 



and 
where 



8 = 1 



/ Ol + 1 O2 + 1 03 04 \ 

10 

Ci c 2 c 3 + 1 c 4 + 1 

\ 1 / 



Spectr(TiT 2 ) = {1,1, * + /,«-/}, 

fli + c 3 j. ^/ a i + c 3 + 4cia 3 — 2aic 3 



/ = 



The straightforward calculation by using (6.7) and (6.8) shows that the eigenvalues of the matrix TiT 2 
are equal to {1,1,1,1}. According to (5.2) these must be the eigenvalues of the matrix T^ which is in 
contradiction to c) of Lemma 5.1. 

Case(2). Assume that even one from the invariants J\, J 2 depends on x\ or x 3 . Let, for example 



(6.9) 



S*°- 



It is useful to consider two additional linear operators 8\ = [5, A] and <5 2 = — [8, 81]. 
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One has 

d d d d d d 

Si = /it h k 1" /37 1" fan — ' * 2 = ( &lX2 + & 32;4)t h (dl^ + ^4)7; — , 

axi 0x2 0x3 0x4 ox\ 0x3 

where 

/g 10 x /1 = -61^1 + (01 - 62)2:2 - b 3 x 3 + (a 3 - 64)2:4, h = ^1^2 + 63214, 

/3 = —d\X\ + (ci - ^2)2:2 - ^32:3 + (03 - ^4)2:4, /i = rfio;2 + d 3 Xi. ' 

We deduce from 5Ji = AJi = that 

<JiJ< = 0, <5 2 J 4 = 0, * = 1,2. 

Consider the partial differential equation 6 J = 0. Solving it one finds that J = K(Yi,Y2,Y 3 ) where 
K{yi 1 y2,y 3 ) is an arbitrary function and 

(6.11) Yi=x 2 , Y 2 = x i7 Y 3 = x A xi - x 2 x 3 . 

d J 

Therefore, in view of (6.9), (6.11) we have J\ — J\(Y]_,Y2,Y$) and — — ^ 0. 

oY 3 

Consequently, as S^Yi = 62Y2 = 0, one gets 

xt d Jl x v _l d Jl x v _l d Jl x v dJl x v 
h.h = w hY + w hY2 + Ws S 2 Y 3 = —S 2 Y 3 . 

This implies 

(6.12) S 2 Y 3 = 0. 

By substituting in (6.12) the expression for Y 3 given by (6.11) we arrive to 

(6.13) b 3 = d x = 0, b 1= d 3 = p, 

for some p G C 

We now use the equation SiJ = which can be written as 

(6.14) Sl J = ^Y + §- Sl Y 2 + §- Sl Y 3 = 0, 

One can show that 

SiY x = P Y l7 

SiY 2 = P Y 2 , 

5iY 3 = v 1 Y 1 2 + V2Y 2 2 +v 3 Y 1 Y2, 

where v\ — c?2 — Ci, v 2 = a 3 — 64, v 3 = 01 — 62 — C3 + (£4. 
Hence, (6.14) yields 

pY m + pY *w2 + ^ + V * Y > + v * Y ^w 3 = °- 

This equation possesses two rational functionally independent solutions Ji(Yi, Y 2 , I3), ^(^i, Y 2 ,Y 3 ) only 
if 

p = 0, «i = u 2 = U3 = 0. 

which gives 

(6.15) a 1 = e 1 +b 2 , c 3 = ei + d 4 , Ci = d 2 = Ci> «3 = &4 = C2, ei,Ci,C2€<C. 
After substitutions of (6.13), (6.15) in (6.2) the matrix R is written as 

/ b 2 + ei a 2 C2 04 \ 

6 2 C2 

Ci c 2 c? 4 + ei c 4 

V Ci d 4 J 
12 



i? = 



Now, consider the characteristic polynomial P(X) of R 

P(X) = A 4 + PxA 3 + P 2 \ 2 + P a X + P 4 , 

where 

Pi = -2(6 2 + di + ei), 

P2 = 36 2 ei - 2C1C2 + 3eirf 4 + 4M 4 + b\ + e\ + d\, 
P3 = ~{d A + b 2 + e 1 )(2b 2 d 4 + b 2 e x + d 4 e 1 - 2(1(2), 
Pi = (b 2 di - C1C2XM4 + fad + d i e 1 - C1C2 + 4)- 

As above, in view of (6.1), (6.2) all eigenvalues of R must be equal to and therefore Pi = 0, 1 < i < 4. 
This system gives 

ei=0, b 2 =m, d,4 = -r)i, r?i+CiC2 = 0, 



and the monodromy matrix T 2 becomes 



T 2 = 



/ 771 + 1 a 2 C2 a 4 \ 

r/i + 1 C2 

Ci c 2 1 - r/i c 4 

V Ci 1 - m J 



The matrix T1T2 has the eigenvalues {1,1,1,1} which contradicts to c) of Lemma 5.1 and proves our 
claim. □ 

Due to our definition of intcgrability we deduce from Theorem 6.2 the following 

Theorem 6.3 The planar three-body problem is meromorphically non-integrable near the Lagrangian parabolic 
solution. 

7. Final remarks 

In the end of 19th century Poincare [18] indicated some qualitative phenomena in the behavior of phase 
trajectories which prevent the appearance of new integrals of a Hamiltonian system besides those which are 
present, but fail to form a set sufficient for complete intcgrability. 

Let M 2n be the phase space, and H : M 2n — > R, H = H + eHi + 0(e 2 ) the Hamiltonian function. 
Suppose that for e = the corresponding Hamiltonian system has an m-dimcnsional hyperbolic invariant 
torus T™. According to the Graff's theorem [5], for small e the perturbed system has an invariant hyperbolic 
torus T™ depending analytically one. It can be shown that TJ™ has asymptotic invariant manifolds A + and 
A~ filled with trajectories which tend to the torus T™ as t — ► +00 and t — ► —00 respectively. In integrable 
Hamiltonian systems such manifolds (called also separatrices), as a rule, coincide. In the nonintegrable cases, 
the situation is different: asymptotic surfaces can have transverse intersection forming a complicated tangle 
which prevent the appearance of new integrals. For a modern presentation of these results see, for example, 

[7]- 

The method of splitting of asymptotic surfaces was applied to the three-body problem by many authors. 

In his book [13] J.K. Moscr described a technique which use the symbolic dynamics associated with a 
transverse homoclinic point. Applying this method, it was shown in [9] that under certain assumptions the 
planar circular restricted three-body problem does not possess an additional real analytic integral. The 
similar result for the Sitnikov problem and the collincar three body problem can be found in [13], [10]. The 
existence and the transverse intersection of stable and unstable manifolds along some periodic orbits in the 
planar three-body problem where two masses are sufficiently small was established in [11], using the results 
obtained in [12]. 

It is necessary to note that Theorem 6.3 implies the nonexistence of a complete set of complex analytic 
first integrals for the general planar three-body problem. To prove the nonexistence of real analytic integrals 
one should use some heteroclinic phenomena and can propose the following line of reasoning : Let M°° be the 
infinity manifold, then the taken Lagrangian parabolic orbit is biasymptotic to it. This is a weakly hyperbolic 
invariant manifold and the reference orbit is a heteroclinic orbit to different periodic orbits sitting in M°°. 
The dynamical interpretation of Theorem 6.3 seems to be the transversality of the invariant stable and 
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unstable manifolds of M°° , along this orbit. A combination of passages near several of these orbits (there 
is all the family obtained by rotation) should allow to prove the existence of a heteroclinic chain. This, in 
turn, gives rise to an embedding of a suitable subshift, with lack of predictability, chaos and implies the 
nonexistence of real analytic integrals. 
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Appendix A. The matrices Aoo, A, B, R, J 

Aqc = Aooij, 1 < i,j < 4. 
_ 1 12a + 5/3 + 5/3a 2 + 26a/3 + 12a 2 3 \/3 (a + l)/?(a - 1) 

Aca, 11 — T n i ^oo, 12 — T 5 : 

c 

A,o,13 = -2 2 3 — , Aoo,14=0, 

i 2 o 3 a 4 p 
3 \/3(a + l)/3(a-l) 1 -12 a + /3 + /3 a 2 - 2 a (3 - 12 a 2 

^4-oc,21 — 7 p; , ^4oo,22 = 



4 aSi 4 aSi 

c 

-4oc,23=0, A OCj2 4 = -2 2 3 — , 

i 2 o 3 a* p 

__ 1 a 2 /? i 3 3 (a + 1) i 2 3 (2 a + 13 /? + 13 (3 a 2 + 24 a /? + 2 a 2 ) 
4»,31 - g 7T3 ■ 

3 V3(/3 + 2a + 4a/3 + /?a 2 + 2a 2 ) (a- 1) /3a 2 i 3 3 i 2 3 



A 



oc.32 



8 i x 3 

1 /3(5a 2 + 14a + 5) 3 \/3 (a + l)/?(a - 1) 

^-oc,33 — ~—. T, ^-oc,34 



^■oo,41 — 

^4oo.42 = 



4 a Si "' 4 a Si 

3 V 



3 V3 (/? + 2a + 4a/3 + /?a 2 + 2a 2 ) (a -1) /3a 2 i 3 3 i 2 3 



8 Si 3 

1 a 2 /3S 3 3 (a + l)S 2 3 (-10a + 7/3 + 7/3a 2 -12a/3-10a 2 ) 

8 S[ 



3 



3 V3(a + l)/3(a-l) 1 /3(10a + a 2 + l) 

Aoo < 43 -~4 a^ ' Aoo ' 44 -I ^ ' 

A = (A tJ ), 1 < *, j < 4. 

1 (a + l)(a/? + 4a + /3) 1 \/3 (a + 1) /? (a - 1) 

An = ~7 e ' ^ 12 — 7 e • 

4 aii 4 aii 

A fr 4 _n 4 _ 1 V3(a + l)/3(a-l) 

-"■13 — ^ t q ; i ^ll4 — U, JT-21 — T T, i 

S 2 2 S 3 3 a A (3 4 aii 

1 10a/3 + 3/3a 2 + 3/3 + 4a 2 +4a ii 

^122 — : , ^124 — 2 



a Si ' ^ i 2 2 S 3 3 a4/3 

. _ 1 a 2 /3 2 i 3 3 (a + l)(a-l) 2 i 2 3 . 1 y/Z (a + l)/3(a - 1) 

^31 — -7 773 , -434 - -7 f; , 

8 Si 4 aii 

, _ n . _ 1 V3(a-l)(a + l) 2 a 2 /3 2 i2 3 S 3 3 _ 1 {a - I) 2 (3 

23 ' 32 8 i! 3 ' 33 4 aii ' 

. 1 V3(a-l)(a + l) 2 a 2 /3 2 i 2 3 i 3 3 . _ 3 a 2 /3 2 i 3 3 (a + l) 3 i 2 3 

-Ail — — 5 , -Ad? — o ■ 

8 ii 3 8 ii 3 

1 V3(a + l)/?(a-l) _3(a + l) 2 /? 

-443 — -7 7, > -444 — 7 5 ■ 

4 aii 4 aii 



1G 



R=(R ij ), l<i,j<4. 



12a + /3 + 6al3 + 2a 2 +l3a 2 _ 1 \/3(a + l)/3(a - 1) 

"11 = _ o C ' "12 = _ o C ■ 

2 aii z aii 

i?13 = 0, i?i4 = 0, 

D 1 \/3(a + l)/3(a-l) 1 (a + 1) (-2a + a/3 + /3) 

-K21 — -r 5 , "22 - x 5 , 

2 aii 2 aii 

"23 = 0, i?24 = 0, 

1 /3a 2 i^ 3 .S3 3 (a + 1) (a 2 + 6 /3a 2 + a + 13a/? + 6/3) 

#31 - -g ^3 



#3 



1 V3 (3 a 2 + 2 /3a 2 + 7a/3 + 3 a + 2/3) (a - l)/3a 2 S 3 3 S 2 3 

S? 



_, _ 1 (a + V3 + 2)(a + 2- v / 3)/3 D _ 1 V3 (a + 1) /3 (a - 1) 

"33 — 7; o ' Hi4 _ o c • 

2 aii z aii 



i?41 — — 
^42 = — 



1 V3 (3 a 2 + 2 /3a 2 + 7a/3 + 3 a + 2/3) (a - 1) /3a 2 i^ 3 i^ 3 
8 S? : 

1 /3a 2 i , 2 3 S , 3 3 (a + l)(-5a 2 + 2/3a 2 -5a-9a/3 + 2/3) 

8 S? ' 



D _ 1 V3(a + l)/3(a-l) 1 (a + \/3 + 2) (a + 2 - y/S) /3 

"43 - 7; q , "44 - -^ 5 ■ 

2 aii z aii 



1 V3(a + l)/?(a-l) 

-Ai - ~o c ' 

2 aii 



Jl2 = ^ 



J={Jij), 1<*,J<4. 
1 (a + l)(-2a + a/3 + /3) 



Jl3 = 0, Jl4 = 

1 2a + /3 + 6a/3 + 2a 2 + /3a 2 

^21 - x ^ , J22 - _ 

2 aii 2 



a Si 
1 \/3(a + l)/3(a-l) 



a Si 



, _, , _,. T 1 V3(a-l)(a+l) 2 a 2 /3 2 52 3 5 3 3 

J23 = U, J 2 4 — U, J31 — — - —3 , 

4 il 

T _ 1 /3 2 a 2 i 2 3 i , 3 3 (a + l)(a 2 + 4a + l) 

J32 ~I 5? ' 

1 \/3(a + l)/3(a-l) _ 1 2a + /? + 6a/3 + 2a 2 + /3a 2 

^33 - ~ c ' "' 34 ~ ~o C ' 

2 aii z aii 



J41 

Ji2 = T 



1 /3 2 a 2 53 d S 3 d (a + 1) (a 2 + 4a + 1) 



4 _ Si 3 

1 V3 (a - 1) (a + l) 2 a 2 /3 2 S 2 3 S 3 3 

4 ii 3 



T _ 1 (a + l)(-2a + a/3 + /3) _ 1 V3 (a + 1) /3 (a - 1) 

J 43 - ~o c ' 44 ~~ ~0 C ' 

2 aii z aii 
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